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A 3 Hc DRIVEN INSTABILITY NEAR THE FULLY CONVECTIVE BOUNDARY 
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ABSTRACT 

We report on the discovery of an instability in low mass stars just above the threshold (~ 0.35 M Q ) 
where they are expected to be fully convective on the main sequence. Non-equilibrium 3 He burning 
creates a convective core, which is separated from a deep convective envelope by a small radiative 
zone. The steady increase in central 3 He causes the core to grow until it touches the surface convection 
zone, which triggers fully convective episodes in what we call the "convective kissing instability" . 
These episodes lower the central abundance and cause the star to return to a state in which is has a 
separate convective core and envelope. These periodic events eventually cease when the 3 He abundance 
throughout the star is sufficiently high that the star is fully convective, and remains so for the rest of its 
main sequence lifetime. The episodes correspond to few percent changes in radius and luminosity, over 
Myr to Gyr timescales. We discuss the physics of the instability, as well as prospects for detecting its 
signatures in open clusters and wide binaries. Secondary stars in cataclysmic variables (CVs) will pass 
through this mass range, and this instability could be related to the observed paucity of such systems 
for periods between two and three hours. We demonstrate that the instability can be generated for 
CV secondaries with mass-loss rates of interest for such systems, and discuss potential implications. 
Subject headings: starsunteriors — stars:evolution — stars: low-mass — binaries: close 



1. INTRODUCTION 

Stars are surprisingly stable for most of their evolution. 
The exceptions provide windows into important transi- 
tions in stellar properties. The helium flash, for example, 
marks the temporal transition from shell hydrogen burn- 
ing to core helium burning. Interiors theory also predicts 
changes in structure as a function of mass. A classic 
example is the prediction that progressively lower mass 
models should have progressively deeper surface convec- 
tion zones, culminating in a transition to a fully convec- 
tive state at around a third of a solar mass. The global 
behavior of models near this transition is traditionally 
thought to be smooth. 

We have uncovered a novel instability for models just 
above the fully convective boundary that has, to the best 
of our knowledge, never before been noted for these ob- 
jects. Non-equilibrium 3 He burning in a marginally sta- 
ble and small radiative core leads to a steady increase 
in 3 He and energy generation, and the development of a 
small convective core separated from the convective enve- 
lope by a radiative buffer zone. Above a critical threshold 
a transition to a fully convective state occurs. Convective 
mixing lowers the core 3 He abundance, quenching the nu- 
clear reactions and causing the core to become radiative 
again. The sudden drop in energy generation leads to a 
sudden drop in radius, followed by a gradual recovery as 
3 He production in the core resumes. This phenomenon 
repeats in a series of pulses, with an amplitude which is 
steadily damped as the envelope 3 He is enriched and the 
drop in core 3 He becomes smaller. For non-interacting 
stars this will apply to a narrow range of masses just 
above the fully convective boundary. The secondary stars 
of cataclysmic variables, however, will all pass through 
this mass range, and the onset of the instability occurs at 
the upper end of the CV period gap. The astrophysical 
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consequences may therefore be surprisingly broad, and 
the instability can be triggered for mass-loss rates of in- 
terest in the CV context. In this paper we describe the 
physics of the instability and present some preliminary 
assessment of its effects. We describe our stellar models 
in Section [21 discuss the details of the physical processes 
responsible for the instability in Section [31 observational 
implications of the instability in Section[4l and close with 
our conclusions in Section [5l 
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2. STELLAR MODELS 

utilize the Yale Rotating Stellar Evolu- 
Code (YREC) (iPinsonneault et al, 
m [199S 



Bahca ll fc Pinsonneault] Il992t iBahcall et al.l 



1989; 



1995 



200 lh to produce a grid of stellar models with a stan- 



dard set of input physics. We create a dense grid near 
the fully convective boundary to investigate the mass 
range over which this instability operates. Models 
range in mass from 0.2 — 0.5 M Q and are spaced every 
0.001 M & . Models generally contain ~ 500 shells and 
take ~ 5000 timesteps reach to 14 Gyr, at which point 
the evolution is truncated. All models include helium 
and heavy elem e nt dif fusion following the procedure 
of iThoul et al.1 (11994D. a t mosph ere and boundary 
conditio ns of | Allard et ah (120 00) , nuclear reaction 
rates of [Adelberger et al.l ( 2011l )~with weak screening 
(|Salpeterlll954fl 1 and employ a mix ing length theory of 
convection ()Cox|[l968t lVitenselll953D. Opacities ar e from 
the Opacity Proje ct (OP) ([Mendoza et al.l 120071) for a 
iGrevesse fc Sauval (|1998[ ) solar mixture , supplemented 
with the low temperature opacities of Ferguson et al.l 
(120051). We utilize the 2006 OPAL equation of state 
(IRoeers et all Il996t iRogers fc Navfonovl |2002[ ) and the 
ISaumon et al.l (|1995t ) EOS for temperature and density 
combinations outside of the OPAL tables. Convec- 
ti ve instability is determined using the Schwarzschild 
ISiterion. We impose an initial 3 He mass fraction of 
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2.95 x 1CT 5 , for a iGrevesse fc Sauvall (|1998|) mixture 
with t he He 3 /He 4 isotopic ratio from I Anders fc Grevessd 
(1989) on a series of pre-main sequence starting models. 
Because pp chain burning is by far the most significant 
source of 3 He, the exact value of this initial abundance 
is of relatively little importance. 

3. 3 He INSTABILITY 

We find that low mass objects at the fully convective 
boundary undergo non-equilibrium 3 He burning during 
the early stages of their main-sequence lifetimes which 
gives rise to periodic fully convective episodes. We dis- 
cuss the nature of the instability, the conditions under 
which it operates, and compare our results obtained with 
YREC to those obtained with the MESA stellar evolu- 
tion code. 

3.1. Nature of the Instability 

3 Hc is produced in low mass stars via proton capture 
on deuterium; the former can be either primordial or 
produced from the pp chain. The pp reaction itself has 
a low cross section but a weak temperature dependence, 
so 3 He can be generated even in relatively cool environ- 
ments given sufficient time. 3 He is destroyed by both 
3 He + 3 He (pp I) and 3 He + 4 He reactions (pp II and 
pp III); the former is far more important at low tem- 
peratures. At high temperatures equilibrium is achieved 
with a 3 He abundance that decreases as the temperature 
increases, and net destruction of 3 He can occur. How- 
ever, at low temperatures the rate of product ion exceeds 
the rate of destruction and 3 He is produced (llbenl|196 
IRood et all H976I: see also [Bocsgaard & SteigmarJ 1198 
for a discussion). One therefore expects a 3 He abundance 
pea k in the outer layers of hi gher mass stars (see Figure 
1 in lPinsonneault et al.l tl989. for an example). For lower 
mass stars the peak shifts towards the center of the star 
and the development of a deep convective envelope ho- 
mogenizes the outer layers. For sufficiently cool central 
temperatures 3 He is never in equilibrium. Therefore the 
nuclear energy generation from 3 He + 3 He burning will 
depend sensitively on the integrated 3 He production over 
the lifetime of the star. 

If there is a radiative core it is possible to have a local 
enhancement in the 3 He abundance which is much larger 
than the average for a fully convective star. This special 
circumstance can therefore create a small convective core, 
even in stars far too cool to support the CNO burning 
responsible for convective cores on the upper main se- 
quence. The growth in central 3 He will cause this core 
to grow in mass. For a sufficiently small radiative core it 
is therefore possible to induce a "convective kissing" in- 
stability where the central and surface convective regions 
merge. 

To understand the nature of the this instability, we 
introduce Figure [TJ which contains plots of the location 
of the base of the envelope convection zone, core 3 He 
fraction (A" 3 jj e ), luminosity, and radius of several stel- 
lar models of different masses. Fully convective episodes 
(when R cz — 0) are in phase with sudden drops in A 3 jj c 
in the core. After an episode the fraction of 3 He in the 
core slowly recovers, along with the luminosity from 3 He 
burning reactions, whose rates are inflated by the pres- 
ence of the elevated 3 He abundances. The total lumi- 
nosity and radius of the star both increase on nuclear 
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Fig. 1. — Behavior of models of several characteristic masses 
through convective episodes. The panels show radius, luminosity, 
location of the base of the convective envelope, and core He 3 frac- 
tion from top panel to bottom panel in order, each as a function 
of time. Curves for YREC models are as follows: solid purple: 
0.35 M model, blue dotted: 0.37 M , green dashed: 0.38 M Q 
(which is too massive to show convective episodes). The red dot- 
dashed curve shows a MESA 0.37 M model. 

timescales in response to the additional burning. The 
response of the core and envelope convection zones is 
shown in Figure [2l which plots the temperature gradi- 
ents of a single 0.360 M Q model at different stages in 
the convective cycles. At times early in the cycle ((a), 
bottom panel of Figure \2§ , the core and envelope con- 
vection zones are separated significantly in radius. As 
time passes and the mean 3 He of the entire star increases 
with each mixing episode, the increased 3 He burning in- 
flates the radiative temperature gradient, resulting in a 
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decrease in the separation between the envelope and core 
convection zones, as seen at times (b) and (c). By time 
(d) the additional luminosity is sufficient to make the 
convection zones meet, at which point the star is fully 
convective, and mixes the 3 He in the core throughout 
the entire model on a convective overturn timescale. In 
the case of a 0.36M Q model, the luminosity from the re- 
actions 3 He + 3 He and 3 He + 4 He in the convective core 
represents between 10 and 35% of the total model lumi- 
nosity, depending on the location within the cycle. The 
dilution of the 3 He in the core results in decreased lumi- 
nosity and thermal support, and the model radius and 
luminosity decrease on a thermal timescale. The sud- 
den loss of substantial 3 He burning is somewhat counter- 
balanced by energy release due to gravitational collapse, 
which results in modest, few percent changes in the total 
luminosity over the course of these episodes. The struc- 
ture returns to a convective core, radiative buffer zone, 
and convective envelope, while the 3 He again begins to 
accumulate in the core. 

The middle panel of Figure [2] shows the temperature 
gradient at the timcstep immediately following the fully 
convective episodes for several cycles. Each successive 
cycle begins with higher 3 He abundances throughout the 
model, since the products of previous cycles are mixed 
throughout the star. Because the total 3 He abundance is 
higher at the onset than in previous cycles, the bound- 
aries of the convective core and envelope begin the cycle 
more closely spaced in radius, and the duration of the 
cycles decreases. Cycles continue until the 3 He abun- 
dance throughout the star is high enough that the model 
is fully convective for the rest of its MS lifetime. 

Figure Q] demonstrates the response of the convection 
zone, central 3 He abundance, radius and luminosity for 
models of different masses over the course of many con- 
vective cycles. More massive models tend to have cy- 
cles of longer duration and greater amplitude, because 
they are further from the fully convective boundary, and 
thus require more 3 He burning to cause them to be con- 
vectively unstable. Likewise, lower mass models have 
shorter cycles with smaller amplitude variations in core 
3 He fraction, luminosity, and radius. In our models (at 
solar composition X, = 0.710342, Z t = 0.018338, with a 
solar calibrated mixing length parameter a = 1.924897), 
the instability occurs over the mass range 0.322 Mq < 
M < 0.365 M Q . For masses 0.266 M Q < M < 0.322 M Q 
the star undergoes a brief period in which it has a ra- 
diative zone between a convective core and convective 
envelope without any episodic mixing, but is then fully 
convective on the remainder of the MS. Models with 
masses M < 0.266 M Q are fully convective (in that they 
never have a radiative core) throughout their MS life- 
times. This fully convective boundary is somewhat lower 
than the typical values quoted in the literature primarily 
due to our definition: we consider the fully convective 
boundary to be the division between models that host a 
radiative core at some point in their lives, but that may 
be fully convective for the vast majority of their lifetimes, 
versus models that never have even a transient radiative 
zone. 

3.2. Checks and Physics 

To verify that this instability is visible in other stel- 
lar evolution codes apart from YREC, we made use 
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Fig. 2. — Temperature gradients and core 3 He abundance in a 
M = 0.360 Mq model. The bottom two panels show the tem- 
perature gradients at different times during the evolution of the 
model, which are labeled on the top panel on a plot of the core 
3 He abundance. In the bottom two panels, the adiabatic temper- 
ature gradient is plotted as a thick dashed black line. Any curves 
that fall below this line represent radiative regions, and any curves 
that follow it convection zones. 



of the publically available M ESA (jPaxton et all I2011L 
|http: / / mesa.sourceforge.net /[ ) stellar evolution code. 
We ran solar calibrated models with Z = 0.016498, 
Y = 0.26275, mi xing length parameter a = 2.01 and 
Z/X = 0.02289 (IGrevesse fc Sauvall IT991 at 4.57 Gyr 
(jBahcall et all 11995ft and the default, out-of-the-box 
physics and numerical parameters released with MESA 
version 3372, with the exception of the model timesteps, 
which we forced to be of order 10 6 years or smaller. We 
found the same qualitative behavior with 3 He at the 
convective boundary in the MESA models, with some 
quantitative differences. The fully convective bound- 
ary in these models occurs for masses M < 0.281 Mq, 
and the mass range over which the convective episodes 
are present is 0.343 M < M < 0.375 M (versus 
0.322 M Q < M < 0.365 M Q for the YREC models). 
We expect that this difference is due to the detailed dif- 
ferences between the input physics in the two different 
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codes. Although the timescales differ somewhat, the be- 
havior is qualitatively the same, which is evidence in fa- 
vor of a physical phenomenon rather than numerical ar- 
tifact. A MESA model for M = 0.37 M is shown in 
Figure Q] among other YREC models for comparison. 

The exact masses at which the star is unstable to 
these convective episodes is also mode l dependent wi thin 
YREC. For example, models run with lKuruczl (|1997f ) at- 
mospheres rather than Allard atmospheres become un- 
stable at M = 0.332 M Q rather than 0.322 M Q . How- 
ever, in both cases the instability occurs over a mass 
range of AM = 0.043 M Q , consistent with a simple 
shift in the location of the fully convective boundary. 
Models run at non-solar metallicities also display un- 
stable behavior over slightly different mass ranges. For 
example, a 0.34 M Q model with [Z/X] = —0.4 (refer- 
enced to solar metallicity) displays instabilities with the 
longer timescales and larger amplitudes of a more mas- 
sive ~ 0.36 solar metallicity object. 

Physical processes such as convective overshoot have 
not been considered here, and will undoubtedly affect 
the exact location in mass of the instability boundaries. 
There are two distinct cases in which overshooting may 
affect the instability; the first is simply by shifting the ef- 
fective location of the boundaries of the convection zones, 
resulting in a shift in the masses at which the instability 
occurs. The second effect is more subtle, and is related 
to the interplay and exchange of material between two 
closely separated convection zones due to overshooting. 

The importance of overshooting in a system with 
closely separated convection zones has been considered 
in the case of A stars, which are thought to have two 
thin surface convection zones, one at the hydrogen ion- 
ization zone, and another at the deeper second ioniza- 
tion of He. These are formally separated by a radiative 
buffer, although it is unclear whether they are truly dis- 
ti nct from one another when overshooting is considered, 
as lFrevtag et all (|1996D ; lLatour et all (|198lD suggest. In 
the particular case we consider here, overshooting would 
serve to connect the core and envelope convection zones 
either chemically or physically before they formally touch 
according to the Schwarzschild criterion. YREC treats 
overshooting as "overmixing" only, and does not adjust 
the temperature gradients in the overshoot regions, and 
so is incapable of capturing the full behavior of the sys- 
tem when the CZs are in this close configuration. One 
could imagine that the dilution of the core 3 He due to 
3 He poor material mixed in through overshooting would 
lead to an overall drop in the 3 He luminosity in the core, 
a global decrease in the radiative temperature gradient, 
and therefore a shrinking of the convective core. Like- 
wise, the mixing of 3 He rich material into the convective 
envelope might lead to increased burning, and cause the 
envelope convection zone to deepen slightly. Essentially, 
if the production timescale for 3 He in the core is much 
shorter than the timescale for the mixing of 3 He poor 
material into the central CZ, then the effects of over- 
shooting should be minimal. If the mixing timescale is 
comparable or shorter than the production timescale for 
3 He, then it is possible that the effects of overshooting 
could prevent these fully convective episodes from occur- 
ring at all. However, if the temperature gradients in the 
overshoot regions are very close to adiabatic, the bound- 
aries of the the two CZs are actually closer in the case 



of overshooting than they would otherwise be. Further 
work on this interaction is therefore of interest. 

We offer, however, that the convective veloci- 
ties in these stars should be low, and the stel- 
lar material strongly stratified. The extent of con- 
vective envelope overshooting in the Sun inferred 
from hclio seismology is less than 0.05 pressure scale 

heights (IBasul 119971). and bo th st udies of cluster 

CMDs dKozhurina-Platais et all 119971; iSaraiedini et al.l 
ll999HVandenBerg fc Stetsonll2004h and asteroseismology 
(jBriquet et al.l 120111 ) suggest that the typical extent of 
core overshooting in intermediate mass stars is 0.1-0.2 
pressure scale heights. There has been very little ob- 
servational or theoretical investigation of the amount of 
overshooting in stars of a few tenths of a solar mass, but 
we might argue, given what information we have about 
other masses, that the degree of overshooting will be min- 
imal. 

Likewise, we have neglected to include the effects of 
magnetic fields, which could be significant for very active 
M-dwarfs, and could potentially stabilize regions of the 
mode l against convection (as in iMullan fc MacDonaldl 
[200lh . 

4. OBSERVATIONAL IMPLICATIONS 
4.1. Single Stars 

To investigate whether the luminosity and radius 
changes associated with the instability are observable, we 
focus on open clusters, in which we have a single, coeval, 
chemically homogeneous sample of stars. We chose a sim- 
ple initial mass function, with dn/dM — AM~ a , where 
A is a normalization constant, and a = 1.3 (|Kroupal 
I2001D . We sample the IMF in 0.001 M Q bins, and deter- 
mine our model luminosities and effective temperatures 
at fixed time for each mass. We add noise to the model 
HR diagram by assigning luminosities and effective tem- 
peratures randomly sampled from a normal distribution 
centered on the theoretical values for each stellar mass. 
We mimic the presence of binaries by randomly selecting 
half of the stars in the theoretical open cluster, and as- 
signing a companion of mass ^-^primary wnere bis a ran- 
dom number between and 1 drawn from a uniform dis- 
tributio n, motivated by the dis tribution of binary mass 
ratios in iRaghavan et al.l (|2010D . We assign the binaries 
the luminosity weighted mean effective temperature of 
the two components, and add the two luminosities. We 
do not account for the effects of activity. The presence of 
the instability produces some weak structure along the 
single-star MS, which will be challenging to detect in any 
open cluster survey. However, a few stars lie below the 
MS, where, with the addition of proper motion or RV 
membership information, they might be recognized as 
'odd' cluster members. Detection of such sources would 
rely on there being an open cluster with enough mem- 
bers, of an interesting age and metallicity such that this 
instability is visible. 

The presence of this instability could also manifest it- 
self as a scatter in radius as a function of age in M dwarf 
wide binary systems in which the stellar parameters are 
well known. Given a sample of stars of the same mass 
and composition with a variety of ages, one should ob- 
serve more scatter in the radii of the younger objects in 
comparison to the older objects. The effect will only be 
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visible if one can disentangle the contributions from stel- 
lar activity and observational uncertainty, both of which 
are current challenges. It is important to note, however, 
that this is potentially an additional source of scatter in 
M-dwarf radii which is not related to activity or environ- 
ment, but rather nuclear burning processes. 

4.2. Cataclysmic variables 

Cataclysmic variables are short period binaries in 
which a white dwarf accretes mass through Roche lobe 
overflow from a relatively unevolved donor star. CVs 
are observed over a wide range of binary periods, rang- 
ing from ~ 8 hrs to just over 1 hour. The primary 
means of angular momentum loss are postulated to be 
magnetic braking associated with the stellar wind of the 
donor star, and losses through gravitational radiation, 
which operates most efficiently for the closest orbital 
separations. One of the most notable features in the 
CV period distrib ution is the so-called " period gap" (see 
iRappaport et all 119831 iSpruit fc Ritterl Tl983) in which 
we observe a paucity of systems with periods of 2-3 hrs. 

This instability is interesting in the context of CV 
systems. Given sufficient time and mass-loss rates, all 
CV donor stars will pass through the mass range over 
which these convective episodes occur in single-star mod- 
els. The response of the donor star radius to mass-loss 
is a critical element of the CV picture; sudden changes 
could truncate the mass-loss by bringing the system out 
of Roche lobe contact, or conversely result in temporarily 
inflated accretion rates. Although a detailed discussion 
of the manner in which this instability interacts with the 
process of regulated mass-loss through Roche-lobe over- 
flow is beyond the scope of this article, we do consider 
the response of single star models to constant mass-loss 
rates. Under the right conditions, this instability induces 
a significant, complicated response in the stellar radius, 
at an extremely interesting mass range in CV systems. 

We apply constant mass-loss rates between 10~ 8 — 
10~ 12 M Q yr _1 to models of 0.4 M Q and a single mass 
loss rate of 5 x 1CP 11 M Q yr _1 to models over a range of 
masses (0.4 — 0.7 M Q ), all of which have been evolved to 1 
Gyr prior to beginning mass-loss. This exploration of of 
mass- loss rates is presented in Figure [31 In general, rela- 
tively unevolved models of the lowest masses (~ 0.4 M©) 
coupled with low mass-loss rates (10 -10 — 10~ n M Q yr _1 ) 
show the largest radius "glitches" due to convective kiss- 
ing. 

Whether or not strong radius glitches appear in these 
models depends sensitively on a combination of the mass- 
loss rate, initial mass, and evolutionary state of the 
model. At fixed initial mass, the instability does not 
manifest itself in models until the mass-loss timescale is 
greater than the thermal timescale, M/M ~ M 2 /RL. 
For high mass loss rates, the mass loss timescale is al- 
ways smaller than the thermal timescale, the star is out 
of thermal balance, and no instability is present, even 
over the mass ranges where we expect it to occur. For 
lower mass-loss rates, the instability is present and the 
glitch in radius becomes visible (see Figure [3]). 

The initial mass of the model sets the 3 He profile and 
abundance, which in turn determines if the 3 He abun- 
dance in the core is sufficient upon arrival at ~ 0.35 M© 
to become unstable. Low mass models have 3 He pro- 
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Fig. 3. — Zoo of mass vs. radius curves for a sampling of mass-loss 
rates obtained using a starting model of 0.4 Mq evolved to 1 Gyr 
before mass-loss is initiated. All mass-loss rates are in solar masses 
per year. The curves correspond to the following mass loss rates: 
dotted (purple): 5 X lCT 10 M yr -1 , dashed (blue): 1 X 10 -10 , 
solid (green): 5 X 10 — , dot-dashed (orange): 1 X 10 — , triple- 
dot-dashed (red): 5 X 10 -12 . Note that the curve for 5 X 10~ 12 is 
truncated at 14 Gyr. 

files peaked closer to the central region, and in general 
the absolute central 3 He abundance increases with de- 
creasing mass for models of the same initial composi- 
tions. Higher mass models, which are more effective at 
burning 3 He, have 3 He profiles that peak at larger radii, 
and generally have more 3 He depleted cores. Finally, be- 
cause models of different masses have different nuclear 
equilibrium timescales for 3 He, in general, the later the 
mass-loss begins, the less 3 He is present in the core of the 
model. These three effects conspire to make the radius 
glitches most visible in relatively unevolved, low mass 
models with low mass-loss rates. We note a small radius 
glitch occurs even in models with low mass-loss rates but 
insufficient 3 He to induce the instability, which occurs as 
3 He is homogenized throughout the star as the model 
crosses the fully convecti ve boundary; we believe this is 
the same feature noted inlD'Antona fc Mazzi telHI (|1982f ) 
and lAndronov fc PinsonneaultJ (|2004f ). 

Apart from the question of whether this instability af- 
fects the degree of mass loss through physically altering 
the radius of the donor star, one should also note that the 
composition of any donated mat erial is influenced by the 
presence of the in stability. Both lShen fc Bildstenl ( 2009') 
and ISharal (|1980f ) suggest that for 3 He mass fractions 
> 10 -3 the ignition of novae may have some sensitivity 
to the 3 He abundance, and that the burning of accreted 
3 He on the primary can lead white dwarfs that are more 
luminous than their average accretion rates suggest. In 
the case of a 0.4 M model with a constant mass loss 
rate of 5 x 10 _11 Moyr _1 , the surface 3 He abundance 
increases by a factor of 2 over the course of the first con- 
vective episode, from ~5x 10~ 4 to 10~ 3 . This increase 
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in the 3 He abundance of accreted material would appear 
at masses above the fully convective boundary, where it 
would otherwise not exist. 

We have negelected to include rotational mixing in our 
models, which could be non-negligible in the case of such 
close binary systems in which the secondary may be rel- 
atively rapidly rotating due to ti dal locking, which could 
help t o drive me ridional flows (Eddingtonl ll925t ISweetl 
f950; Zah^EMl) and result in mixing. Furthermore, 
mixing may also be driven by tidal distortions of the 
potential due to its close companion (iTassoul fe Tassoull 
Il982t lAndronov fe Pinsonneaultll2004) ). If the timescale 
for mixing is shorter than the timescale in between con- 
vective kissing episodes, the mixing of material could 
serve to eliminate episodic mixing events. With the in- 
clusion of tidal deformation, the estimates of the min- 
imum timescale for mixing are of order the timescale 
between convective episodes for systems at 3 hrs (see 
lAndronov fc Pinsonneault! 12004 Fig. 9 for timescales). 
The initial feature in the curve in Figure [3J as the 3 He 
in the core is first mixed throughout the envelope would 
remain even in the case of strong mixing, but the subse- 
quent variations in the radius would be damped. There- 
fore, the precise shape of these curves may be modified 
in the case in which rotational mixing is important, and 
our neglect of mixing is an important caveat to this dis- 
cussion. 



5. CONCLUSION 

We have discovered a novel instability in low mass, 
M ~ 0.35 M Q , stellar models in which a buildup of 3 He 
in the core results in brief convective episodes within the 
first few Gyrs on the main sequence. These events are 
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accompanied by changes in the radius and luminosity of 
the star. We have verified that the same instability ap- 
pears to occur in the MESA code as well, which suggests 
that it is not a merely numerical artifact. Stars under- 
going these convective episodes may be visible as both 
single stars or as members of binary systems. We have 
demonstrated that the radius of the donor star in a CV 
system could have a sizable response to the onset of the 
convective kissing instability, given donor stars in which 
the mass-loss rates, initial stellar mass, and evolution- 
ary state ensures that there is sufficient 3 He available in 
the core by the time the stellar mass reaches and inter- 
esting range. The complicated structure of the glitches, 
which includes both increases and decreases in radius 
over the course of the instability, could interact in inter- 
esting ways with the process of CV mass-loss, either by 
severing contact, or by leading to temporarily inflated 
mass-loss rates, provided they are not suppressed by ro- 
tational mixing processes. We suggest that this instabil- 
ity be investigated in the context of full CV evolution- 
ary models, where the interaction between the mass-loss 
rate and the convective kissing instability could be better 
quantified. 
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